Los Alamos National Laboratory Poisson Superfish
Program Fish written by
Ron F. Holsinger in collaboration with Klaus Halbach

The original Poisson Superfish codes were developed
by Ron F. Holsinger in collaboration with Klaus Halbach.
These programs are provided as a service to the accelerator
community by the Los Alamos Accelerator Code Group (LAACG).

(c) Copyright 1985-2004, by the Regents of the University of California.        
This software was produced under U. S. Government contract W-7405-ENG-36        
by Los Alamos National Laboratory, which is operated by the University          
of California for the U. S. Department of Energy.  Neither the Government       
nor the University makes any warranty, express or implied, or assumes           
any liability or responsibility for its use, or represents that use of          
this software would not infringe privately owned rights.                        
Unpublished - rights reserved under Copyright Laws of the United States.        


Program Fish  7.12  released 3-25-2004 

Starting from file C:\temp\temp2\My code\Tmp\CEL01.T35.
Program file: C:\CPP\CODEBLOCKS\GENDTL\BIN\RELEASE_QT5\SF_FILES\MDTFISH.EXE
SF.INI file: C:\temp\temp2\My code\Tmp\SF.INI  1-10-2024  13:10:18
Memory used for CEL01.T35 arrays:       2.105 M
Memory used for other Fish arrays:       3.569 M

Problem description:
Calcul_156 DTL(1)
cel0.txt CEL01   1-gap problem
Problem file: C:\temp\temp2\My code\Tmp\CEL01.AM  1-10-2024  13:19:26
Problem file length:        2973 bytes
Originating program: MDTfish    Problem type:  0-DT DTL

Coordinates and lengths have dimensions of centimeters.

Problem constants and variables.
Letter A in the code column indicates a value supplied in the
Automesh input file.

Variable Code         Value     Description
ALPHAT          3.930000000E-03 Temperature coefficient of resistance
ASCALE           3767.30313     Scaling factor for H at drive point
BETA                    0.0     Particle velocity
BETA1           0.100000000     Starting BETA in transit-time table
BETA2           0.950000000     Ending BETA in transit-time table
BETAD                   0.0     Design beta for a cavity
CCLDELK          1.00000000     Increment for coupling for table in SFO
CCLMAXK          6.00000000     Highest coupling for table in SFO
CCLMINK          1.00000000     Lowest coupling for table in SFO
CLENGTH                 0.0     Cavity length for normalization in SFO
CLIGHT          2.997924580E+10 Exact speed of light in cm/sec
CONV             1.00000000     Length conversion (number of units per cm)
DBETA           5.000000000E-02 BETA increment in transit-time table
DELFR                   0.0     Frequency step size for a resonance search
DIAGDLL                   0     If 1, DLL writes diagnostics to DiagDLL.txt
DPHI             180.000000     Phase length of the problem geometry
DSLOPE     A    -1.00000000     Slope of D(k^2) function
DSTOLER         2.000000000E-02 Tolerance required on D(k^2)
DX1             0.100176820     First X mesh interval (at XMIN)
DXMIN           2.500000000E-02 Minimum X mesh interval (found by Automesh)
DYMIN           1.733633442E-02 Minimum Y mesh interval (found by Automesh)
ENORM            1000000.00     Field normalization for NORM=4 option
EPS0            8.854187818E-12 Permittivity of free space
EPSIK      A    1.000000000E-08 Frequency convergence parameter
EPSO       A    1.000000000E-06 Convergence parameter in mesh optimization
EZERO      A     1301967.90     E0 for normalization in SFO when NORM=0
EZEROT           1000000.00     E0*T for normalization in SFO when NORM=1
FMU0            1.256637061E-06 Permeability of free space
FREQ       A     348.985634     RF cavity resonant frequency
FREQD                   0.0     Design frequency for a cavity (MHz)
HPHI             5000.00000     Normalization magnetic field for NORM=2
IBETA                     0     If >0, SFO writes transit-time vs BETA
ICCC                      1     1 for real arrays, 2 for complex arrays
ICCP                      1     If 1, compute material power loss
ICORNER1                  0     First corner segment for computing average H
ICORNER2                  0     Last corner segment for computing average H
ICYLIN     A              1     0 for X,Y problems, 1 for Z,R problems
IMAX                    136     KMAX+2
INFODATA                 14     Number of tuning-code parameters
IOBSEG                   -1     First segment of the CCL outer boundary
IPIVOT                    1     Pivoting in matrix inversion routines
IRESID                    0     If 1, calculate potential residuals
IRMAX                    25     Used in optimization of RHOXY
IRTYPE     A              0     Surface resistance option indicator
ISLOT                     0     If 1, SFO computes coupling-slot power loss
ITFILE                    0     If 1, SFO writes transit-time plot file
ITOT                  58480     (KMAX+2)*(LMAX+2)
KDRI                    134     K coordinate of the drive point
KMAX                    134     Number of horizontal logical mesh points
KMETHOD                   0     Wavenumber computation method
KPROB      A              1     Problem type indicator (Superfish)
LAST35                    1     Code for last program to update T35 file
LCYCLE                   92     Iteration number in mesh optimization
LDRI                    279     L coordinate of the drive point
LINT                      1     Logical-mesh coordinate for Ez integration
LMAX                    428     Number of vertical logical mesh points
MAXCY                    19     Maximum number of cycles (-1: use default)
MAXPPR                 1264     Maximum points per region
NAIR                  37853     Number of air points
NBND                    134     Number of Dirichlet boundary points
NBSLF      A              1     Left-side boundary condition
NBSLO      A              0     Lower boundary condition
NBSRT      A              1     Right-side boundary condition
NBSUP      A              1     Upper boundary condition
NDRI                  38079     Drive point index = IRLAX(NPINP)
NEGAT                     0     Zero-area triangle indicator
NFE                       0     Number of iron points
NHSTEM                    1     Number of half stems
NINTER                    0     Number of interface points
NMATR                     0     Number of material records in T35 file
NORM       A              0     Normalization method in SFO
NPBOUND                1263     Total number of boundary points in the mesh
NPINP                 37988     Total points in problem
NPONTS                33373     Number of unknown relaxation points
NREG                     15     Number of regions
NRMSEG                    1     Normalization segment number for NORM=2
NSEG                     54     Number of boundary segments
NSPL                      1     Number of special-potential points
NSTEP                     0     Number of steps for a resonance search
OMEGAM          1.000000000E-03 Used in optimization of RHOXY
PI               3.14159265     The number pi to machine precision
PLCELL     A     360.000000     Phase length per cell for multicell problems
RESIDR          1.000000000E-08 Residual resistance of a superconductor
RFMU             1.00000000     Permeability for rf surface resistance
RHO             1.724100000E-06 Material resistivity (Ohm-m)
RHOR            1.724100000E-06 Reference resistivity (Ohm-cm) at TEMPR
RHOXY            1.60000000     Over-relaxation factor in mesh optimization
RMASS      A    -2.00000000     Rest mass energy of particle in SFO
RS                      0.0     RF surface resistance (Ohms)
RSTEM            1.00000000     Stem radius in cm
SLOSS           3.000000000E-02 Coupling-slot power factor per % coupling
TC               9.20000000     Critical temperature of a superconductor
TEMPC            20.0000000     Normal conductor operating temperature
TEMPK            2.00000000     Operating temperature of a superconductor
TEMPR            20.0000000     Reference temperature for IRTYPE=3
TRIAVG          3.318874850E-03 Average area of all triangles
TRIMAX          9.475319874E-03 Area of the largest positive-area triangle
TRIMIN          2.660097982E-05 Area of the smallest positive-area triangle
VOLUME           23717.4833     Cavity volume (cylindrical symmetry only)
XDRI       A     10.5564650     X coordinate of the drive point
XMAXG            10.5564650     Upper X bound of the problem geometry
XMING                   0.0     Lower X bound of the problem geometry
XNORM1                  0.0     Starting X for NORM=4 integration path
XNORM2                  0.0     Ending X for NORM=4 integration path
XYAREA           247.959778     Total cross sectional area
YDRI       A     6.60606602     Y coordinate of the drive point
YMAXG            27.2200000     Upper Y bound of the problem geometry
YMING                   0.0     Lower Y bound of the problem geometry
YNORM1                  0.0     Starting Y for NORM=4 integration path
YNORM2                  0.0     Ending Y for NORM=4 integration path
ZCTR                    0.0     Reference Z in transit-time integrals
-------------------------------------------------------------------------------
No materials other than air or empty space.
-------------------------------------------------------------------------------
      Kmax =  134, Lmax =  428                       k^2 = 0.005349749
      ----------------------------------------------------------------
      Writing temporary data arrays to memory.
      Memory used:        60.347 M
      ----------------------------------------------------------------
      Total number of mesh points = 58480             freq = 348.98563
                                                 D(k^2) = -1.47406E-07
      Solution time =    0.906 sec.                       End cycle  0
      ----------------------------------------------------------------
                                                        Start cycle  1


               Using a one-point line with slope=  -1.0.


                              del(k^2) = -1.47406E-07  k^2 = 0.0053496
                                                      freq = 348.98083
                                                 D(k^2) = -1.92639E-11
      Solution time =    0.884 sec.                       End cycle  1
      ----------------------------------------------------------------
                                                        Start cycle  2

                         k^2        D(k^2)      D'(k^2)     D"(k^2)
      Current point  0.0053496  -1.9264E-11  -9.9987E-01
      Upper bound    0.0053496  -1.9264E-11  -9.9987E-01

      BEST=1             k^2        D(k^2)      D'(k^2)
      1 (cycle  1)   0.0053496  -1.9264E-11  -9.9987E-01
      2 (cycle  0)   0.0534975  -1.4741E-07  -9.9987E-01


               Using the two-point parabola with slope=-1 at D(k^2)=0.


                              del(k^2) = -1.92639E-11  k^2 = 0.0053496
                                                      freq = 348.98083
      ----------------------------------------------------------------
      Solution converged at the start of cycle 2.
      Initial frequency = 348.98563        Final frequency = 348.98083
      Elapsed time  =    1.85 sec

Problem variables computed by this code.

Variable Code         Value     Description
DKSQ           -1.926391251E-11 Change in k^2 after an iteration in Fish
DSLOPE         -0.999869314     Slope of D(k^2) function
ERG              1750.60073     Integral (H^2 r dr dz)
FREQ             348.980825     RF cavity resonant frequency
ICYCLE                    2     Present iteration number
METHOD                    3     Method used to get frequency in root finder
Q2I                     0.0     1/2Q passed from CFish to SFO
RESIK           3.600999262E-09 Residual = DKSQ/XKSQ
XK0             7.314097256E-02 The wave number k = 2pif/c
XKSQ            5.349601868E-03 Square of the wave number
